eazy-gnuplot 'gnuplot Cookbook' Companion¶
This document is a companion for eazy-gnuplot users going through the 'excellent book by Lee Philips 'gnuplot Cookbook' from Packt Publishing. I highly recommend buying this book.
In this companion we demonstrate implementing the plots in the book using eazy-gnuplot. Liberties have been taken plots may not follow exactly what is shown in the book (when they don't they follow the spirit).
In [1]:
(quicklisp-client:quickload '(:eazy-gnuplot :clml.statistics :clml.utility))
(use-package :eazy-gnuplot)
Out[1]:
1 Plotting Curves Boxes Points and More¶
In [2]:
(defun function-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output)
(func-plot "[-5:5] (sin(1/x) - cos(x))*erfc(x)"))
output)
(png-from-file (function-plot "images/function-plot.png"))
Out[2]:
In [3]:
(defun seperate-y-axis-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output :y2tics '(-100 10) :ytics
:nomirror)
(func-plot "sin(1/x) axis x1y1,100*cos(x) axis x1y2"))
output)
(png-from-file (seperate-y-axis-plot "images/seperate-y-axis-plot.png"))
Out[3]:
In [4]:
(defun seperate-y-axis2-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output :x2tics '(-20 2) :xtics
:nomirror :xrange '("[-10:10]") :x2range '("[-20:10]") :samples
100)
(func-plot "sin(1/x) axis x1y1, 100*cos(x-1) axis x2y2"))
output)
(png-from-file (seperate-y-axis2-plot "images/seperate-y-axis2-plot.png"))
Out[4]:
In [5]:
(defun scatter-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output)
(plot
(lambda ()
(loop for p in (map 'list (lambda (x y) (list x y))
(clml.statistics:rand-n
(clml.statistics:chi-square-distribution 100) 300)
(clml.statistics:rand-n
(clml.statistics:chi-square-distribution 10) 300))
do (format t "~&~{~a~^ ~}" p)))
:with '(:points :pt 7)))
output)
(png-from-file (scatter-plot "images/scatter-plot.png"))
Out[5]:
In [6]:
(defun plotting-boxes-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output :style '(:fill :pattern))
(func-plot "[-6:6] besj0(x) with boxes, sin(x) with boxes"))
output)
(png-from-file (plotting-boxes-plot "images/plotting-boxes-plot.png"))
Out[6]:
In [7]:
(defun plotting-boxes-solid-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output :style '(:fill :solid))
(func-plot "[-6:6] besj0(x) with boxes, sin(x) with boxes"))
output)
(png-from-file
(plotting-boxes-solid-plot "images/plotting-boxes-solid-plot.png"))
Out[7]:
In [8]:
(defun plotting-circles-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output)
(plot
(lambda ()
(loop for p in (map 'list (lambda (x y z)
(declare (ignorable z))
(list x y))
(clml.statistics:rand-n
(clml.statistics:chi-square-distribution 100) 300)
(clml.statistics:rand-n
(clml.statistics:chi-square-distribution 10) 300)
(clml.statistics:rand-n
(clml.statistics:chi-square-distribution 1) 300))
do (format t "~&~{~a~^ ~}" p)))
:with '(:circles)))
output)
(png-from-file (plotting-circles-plot "images/plotting-circles-plot.png"))
Out[8]:
In [9]:
(defun drawing-filled-curves-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output)
(plot
(lambda ()
(loop for i from 0 upto 50
do (format t "~&~a ~a" i (sin i))))
:with '(:filledcurves :above :y1 = 0.07)))
output)
(png-from-file
(drawing-filled-curves-plot "images/drawing-filled-curves-plot.png"))
Out[9]:
In [10]:
(defun drawing-blue-filled-curves (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output :style '(fill pattern 5))
(plot
(lambda ()
(loop for i from 0 upto 50
do (format t "~&~a ~a" i (sin i))))
:lt '(rgb "blue") :with '(:filledcurves :above :y1 = 0.07)))
output)
(png-from-file
(drawing-blue-filled-curves "images/drawing-blue-filled-curves.png"))
Out[10]:
This exercise requires Gnuplot 5
In [11]:
(defun filled-curves-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(:pngcairo) :output output :style
'(fill pattern 5))
(plot
(lambda ()
(loop for i from 0 upto 1 by 0.1
do (format t "~&~a ~a ~a" i i (sin i))))
:using '(1 2 3) :lw (list 3 'notitle) :using '(1 2)
:lw (list 3 'notitle) :using '(1 3) :lw '(3 notitle) :with
'(:filledcurves :above :y1 = 0.07)))
output)
(png-from-file (filled-curves-plot "images/filled-curves-plot.png"))
Out[11]:
In [12]:
(defun finance-bars-plot (output)
(let ((s
(clml.utility.data:fetch
"http://www.google.com/finance/getprices?i=60&p=10d&f=d,o,h,l,c,v&df=cpct&q=ibm"
:stream t)))
(clml.utility.data:process-finance-header s)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output :bars 2)
(plot
(lambda ()
(loop for row across (subseq
(clml.utility.csv:read-csv-stream s :type-spec
'(float float
float
float
float
float))
0 100)
do (format t "~&~{~a~^ ~}" (coerce row 'list))))
:using '(0 2 3 4 5) :with 'financebars)))
output)
(png-from-file (finance-bars-plot "images/finance-bars-plot.png"))
Out[12]:
In [13]:
(defun histogram-like-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output)
(plot
(lambda ()
(loop for i from -2 upto 2 by 0.1
do (format t "~&~a ~a" i (- (* i i)))))
:with '(:histeps)))
output)
(png-from-file (histogram-like-plot "images/histogram-like-plot.png"))
Out[13]:
In [14]:
(defun histogram-stacked-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output :style
'(fill solid 1.0 border lt -1) :style '(data histograms) :style
'(histogram rowstacked))
(plot
(lambda ()
(loop for i from -2 upto 2 by 0.1
do (format t "~&~a ~a ~a ~a" i (- (* i i)) (* i i) i)))
:using (list "(-$2)") :using (list "(20*$2)") :notitle nil))
output)
(png-from-file (histogram-stacked-plot "images/histogram-stacked-plot.png"))
Out[14]:
In [15]:
(defun multiple-histograms-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title
"multiple histograms" :key '(invert reverse left outside) :key
'(autotitle columnheader) :style '(data histogram) :style
'(histogram rowstacked) :style '(fill solid border -1))
(plot
(lambda ()
(loop for r in '((8.01 1 5 1) (8.02 3 5 1) (8.03 4 4 1) (8.04 3 4 1)
(8.05 1 2 1))
do (format t "~&~{~^~a ~}" r)))
:using '(2 "xtic(1)") :title "col0" :using 2 :title "col1" :using 3 :title
"col2" :using 4 :title "col3"))
output)
(png-from-file
(multiple-histograms-plot "images/multiple-histograms-plot.png"))
Out[15]:
In [16]:
(defun error-bars-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title "error bars" :bars 3
:pointsize 3 :xrange (list "[ 8.0:8.06]"))
(plot
(lambda ()
(loop for r in '((8.01 1 0.2) (8.02 3 0.3) (8.03 4 0.4) (8.04 3 0.4)
(8.05 1 0.1))
do (format t "~&~{~^~a ~}" r)))
:using "1:(-$2):3 with errorbars" :title "error" :using "1:(-$2):3 pt 7"
:title "col0"))
output)
(png-from-file (error-bars-plot "images/error-bars-plot.png"))
Out[16]:
In [17]:
(defun boxerror-bars-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title "error bars" :style
'(fill pattern 2 border) :xrange (list "[8.0:8.06]"))
(plot
(lambda ()
(loop for r in '((8.01 1 0.2) (8.02 3 0.3) (8.03 4 0.4) (8.04 3 0.4)
(8.05 1 0.1))
do (format t "~&~{~^~a ~}" r)))
:using "1:($2):3 with boxerrorbars" :title "boxers"))
output)
(png-from-file (boxerror-bars-plot "images/boxerror-bars-plot.png"))
Out[17]:
In [18]:
(defun candlestick-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title "candlestick"
:xrange (list "[0:5]") :yrange (list "[0:25]") :boxwidth
'(0.2 absolute))
(plot
(lambda ()
(loop for r in '((1 5 7 10 15 24 0.3 "'quick'")
(2 6 8 11 16 23 0.4 "'fox'")
(3 5 7 11 17 22 0.5 "'lazy'")
(4 6 9 10 18 21 0.3 "'dog'"))
do (format t "~&~{~^~a ~}" r)))
:using
(list
"1:3:2:6:5:xticlabels(8) with candlesticks title 'quartiles' whiskerbars")
:using (list "1:4:4:4:4 with candlesticks lt -1 notitle")))
output)
(png-from-file (candlestick-plot "images/candlestick-plot.png"))
Out[18]:
In [19]:
(defun impluse-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title "candlestick"
:xrange (list "[0:5]"))
(plot
(lambda ()
(loop for r in '((1 5 7 10 15 24 0.3 "'quick'")
(2 6 8 11 16 23 0.4 "'fox'")
(3 5 7 11 17 22 0.5 "'lazy'")
(4 6 9 10 18 21 0.3 "'dog'"))
do (format t "~&~{~^~a ~}" r)))
:using (list "1:6:xticlabels(8) with impulses")))
output)
(png-from-file (impluse-plot "images/impluse-plot.png"))
Out[19]:
In [20]:
(defun stem-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title "candlestick"
:xrange (list "[0:5]") :yrange (list "[0:25]"))
(plot
(lambda ()
(loop for r in '((1 5 7 10 15 24 0.3 "'quick'")
(2 6 8 11 16 23 0.4 "'fox'")
(3 5 7 11 17 22 0.5 "'lazy'")
(4 6 9 10 18 21 0.3 "'dog'"))
do (format t "~&~{~^~a ~}" r)))
:using '(1 6 "xticlabels(8)") :with '(impulses lw 4 title "counts") :using
'(1 6) :with '(points pt 9 notitle)))
output)
(png-from-file (stem-plot "images/stem-plot.png"))
Out[20]:
In [21]:
(defun parametric-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title "parametric"
:samples 1000 :parametric 'nil :trange (list "[-7:20]"))
(func-plot "sin(7*t), sin(111*t)" :notitle 'nil))
output)
(png-from-file (parametric-plot "images/parametric-plot.png"))
Out[21]:
In [22]:
(defun parametric2-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title "parametric"
:samples 1000 :parametric 'nil :trange (list "[-7:20]"))
(func-plot "sin(7*t), sin(111*t)" :notitle 'nil))
output)
(png-from-file (parametric2-plot "images/parametric2-plot.png"))
Out[22]:
In [23]:
(defun polar-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title "polar" :xtics
'(axis nomirror) :ytics '(axis nomirror) :zeroaxis 'nil :samples
500 :polar 'nil :trange (list "[0:12*pi]"))
(gp :unset :border)
(gp :unset :raxis)
(func-plot "cos(7*t), sin(111*t)"))
output)
(png-from-file (polar-plot "images/polar-plot.png"))
Out[23]:
2 Annotating With Labels and Legends¶
In [24]:
(defun labeling-the-axis-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title
"axis labels and lagends" :parametric 'nil :xlabel "time (sec.)"
:ylabel "amplitude" :y2label "error magnitude")
(plot
(lambda ()
(loop for x from -5 upto 5 by 0.1
do (format t "~&~a ~a ~a" x (sin x) (log x))))
:using '(1 2 "(sgn($2))") :with 'filledcurves :using '(1 2) :with
'(lines notitle) :using '(1 "(sgn($2))") :with '(lines notitle) :using
'(1 "(abs(sgn($2)-$2))") :with '(lines axis x1y2)))
output)
(png-from-file (labeling-the-axis-plot "images/labeling-the-axis-plot.png"))
Out[24]:
This plot generates a postscript file¶
If you have a postscript viewer like ghostscript you open the file at the path given by the output of `setting-label-size-plot. Other wise we have included a conversion of it to a png below
In [25]:
(defun setting-label-size-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(postscript landscape "courier, 18") :output output
:title "axis labels and lagends" :parametric 'nil :xlabel
"time (sec.)" :ylabel "amplitude" :y2label "error magnitude")
(plot
(lambda ()
(loop for x from -5 upto 5 by 0.1
do (format t "~&~a ~a ~a" x (sin x) (log x))))
:using '(1 2 "(sgn($2))") :with 'filledcurves :using '(1 2) :with
'(lines notitle) :using '(1 "(sgn($2))") :with '(lines notitle) :using
'(1 "(abs(sgn($2)-$2))") :with '(lines axis x1y2)))
output)
(setting-label-size-plot "images/setting-label-size-plot.ps")
Out[25]:
In [26]:
(defun adding-a-legend-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title
"axis labels and lagends" :parametric 'nil :xlabel "time (sec.)"
:ylabel "amplitude" :y2label "error magnitude")
(plot
(lambda ()
(loop for x from -5 upto 5 by 0.1
do (format t "~&~a ~a ~a" x (sin x) (log x))))
:using '(1 2 "(sgn($2))") :with '(filledcurves notitle) :using '(1 2)
:with '(lines title "square") :using '(1 "(sgn($2))") :with
'(lines title "fouier") :using '(1 "(abs(sgn($2)-$2))") :with
'(lines axis x1y2 title "error")))
output)
(png-from-file (adding-a-legend-plot "images/adding-a-legend-plot.png"))
Out[26]:
In [27]:
(defun putting-a-box-around-a-legend-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title
"axis labels and lagends" :parametric 'nil :xlabel "time (sec.)"
:ylabel "amplitude" :y2label "error magnitude" :key
'(box lt -1 lw 2 opaque))
(plot
(lambda ()
(loop for x from -5 upto 5 by 0.1
do (format t "~&~a ~a ~a" x (sin x) (log x))))
:using '(1 2 "(sgn($2))") :with '(filledcurves notitle) :using '(1 2)
:with '(lines title "square") :using '(1 "(sgn($2))") :with
'(lines title "fouier") :using '(1 "(abs(sgn($2)-$2))") :with
'(lines axis x1y2 title "error"))
)
output)
(png-from-file
(putting-a-box-around-a-legend-plot
"images/putting-a-box-around-a-legend-plot.png"))
Out[27]:
In [28]:
(defun adding-a-label-with-arrow-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title
"axis labels and lagends" :parametric 'nil :xlabel "time (sec.)"
:ylabel "amplitude" :y2label "error magnitude" :label
'("\"approximation error\"" right at -0.5 #\, 0.9 offset -0.5 #\,
0)
:arrow
'(1 from first -0.5 #\, 0.9 to 0 #\, 1 lt 1 lw 2 front size 0.3
#\, 15))
(plot
(lambda ()
(loop for x from -5 upto 5 by 0.1
do (format t "~&~a ~a ~a" x (sin x) (log x))))
:using '(1 2 "(sgn($2))") :with '(filledcurves notitle) :using '(1 2)
:with '(lines title "square") :using '(1 "(sgn($2))") :with
'(lines title "fouier") :using '(1 "(abs(sgn($2)-$2))") :with
'(lines axis x1y2 title "error")))
output)
(png-from-file
(adding-a-label-with-arrow-plot "images/adding-a-label-with-arrow-plot.png"))
Out[28]:
In [29]:
(defun using-unicode-characters-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :set :samples 500)
(gp :set :encoding 'utf8)
(gp :set :title "favorite graph of ǫrnólfr þórðr")
'font
"helvetica, 24"
(func-plot "[0:10] sin(1/x)"))
output)
(png-from-file
(using-unicode-characters-plot "images/using-unicode-characters-plot.png"))
Out[29]:
3 Applying Colors and Styles¶
In [30]:
(defun coloring-curves-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo solid lw 4))
(gp :set :border '(lw 0.25))
(gp :set :key '(top left))
(func-plot
"[0:1] x**0.5 lc rgb 'orange', x lc rgb 'steelblue',x**2 lc rgb 'bisque', x**3 lc rgb 'seagreen'"))
output)
(png-from-file (coloring-curves-plot "images/coloring-curves-plot.png"))
Out[30]:
In [31]:
(defun styling-curves-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo solid lw 4))
(gp :set :key '(top left))
(gp :set :style '(:line 1 lt 3 lw 4 lc rgb ("'#990042'")))
(gp :set :style '(:line 2 lt 7 lw 3 lc rgb ("'#31f120'")))
(gp :set :style '(:line 3 lt 2 lw 3 lc rgb ("'#0044a5'")))
(gp :set :style '(:line 4 lt 9 lw 4 lc rgb ("'#888888'")))
(func-plot "[0:1] x**0.5 ls 1, x ls 2, x**2 ls 3, x**3 ls 4"))
output)
(png-from-file (styling-curves-plot "images/styling-curves-plot.png"))
Out[31]:
In [32]:
(defun styling-curves--userstyles-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo solid lw 4))
(gp :set :key '(top left))
(gp :set :style '(increment userstyles))
(gp :set :style '(:line 1 lt 3 lw 4 lc rgb ("'#990042'")))
(gp :set :style '(:line 2 lt 7 lw 3 lc rgb ("'#31f120'")))
(gp :set :style '(:line 3 lt 2 lw 3 lc rgb ("'#0044a5'")))
:style
'(:line 4 lt 9 lw 4 lc rgb ("'#888888'"))
(func-plot "[0:1] x**0.5, x, x**2, x**3 "))
output)
(png-from-file
(styling-curves--userstyles-plot
"images/styling-curves--userstyles-plot.png"))
Out[32]:
In [33]:
(defun transparency-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :set :style '(fill transparent solid 0.3))
(func-plot
"[0:pi] sin(x)**2 with filledcurves above y1=0 lc rgb '#00ffff', 0.75*cos(2*x)**2 with filledcurves above y1=0 lc rgb '#aa00aa'"))
output)
(png-from-file (transparency-plot "images/transparency-plot.png"))
Out[33]:
In [34]:
(defun plotting-points-with-curves-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(func-plot "sin(x)/x with linespoints ps 2 pt 6"))
output)
(png-from-file
(plotting-points-with-curves-plot
"images/plotting-points-with-curves-plot.png"))
Out[34]:
In [35]:
(defun changing-the-point-style-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :set :key '(top left))
(gp :set :style '#'linespoints)
(gp :set :style '(line 1 lw 4 lc rgb ("'#990042'") ps 2 pt 6 pi 5))
(gp :set :style '(line 2 lw 3 lc rgb ("'#31f120'") ps 2 pt 12 pi 3))
(gp :set :style '(line 3 lw 3 lc rgb ("'#0044a5'") ps 2 pt 9 pi 5))
(gp :set :style '(line 4 lw 4 lc rgb ("'#888888'") ps 2 pt 7 pi 4))
(func-plot "[0:1] x**0.5 ls 1, x ls 2, x**2 ls 3, x**3 ls 4"))
output)
(png-from-file
(changing-the-point-style-plot "images/changing-the-point-style-plot.png"))
Out[35]:
In [36]:
(defun plotting-with-objects-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo) :object
'(1 circle at graph 0.5 #\, 0.5 size graph 0.25 fillcolor rgb
("'blue'") fillstyle solid arc ("[0:30]") front)
:object
'(2 circle at graph 0.5 #\, 0.5 size graph 0.25 fillcolor rgb
("'red'") fillstyle solid arc ("[30:80]") front)
:object
'(3 circle at graph 0.5 #\, 0.5 size graph 0.25 fillcolor rgb
("'orange'") fillstyle solid arc ("[80:180]") front)
:object
'(4 circle at graph 0.5 #\, 0.5 size graph 0.25 fillcolor rgb
("'green'") fillstyle solid arc ("[180:240]") front)
:object
'(5 circle at graph 0.5 #\, 0.5 size graph 0.25 fillcolor rgb
("'sandybrown'") fillstyle solid arc ("[240:360]") front))
(gp :unset :key)
(gp :unset :tics)
(func-plot "[0:1][0:1] -1"))
output)
(png-from-file
(plotting-with-objects-plot "images/plotting-with-objects-plot.png"))
Out[36]:
4 Controlling Your Tics¶
In [37]:
(defun adding-minor-tics-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :set :grid 'nil)
(gp :set :mxtics 4)
(gp :set :mytics 2)
(func-plot "[0:2*pi] sin(12*x)*exp(-x/4)"))
output)
(png-from-file (adding-minor-tics-plot "images/adding-minor-tics-plot.png"))
Out[37]:
In [38]:
(defun placing-tics-on-second-y-axis-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output
:terminal '(pngcairo)) ;
(gp :set :ytics 'nomirror)
(gp :set :y2tics 0.4)
(gp :set :xtics '("pi/4."))
(func-plot "[0:2*pi] sin(x) axis x1y1, 2*cos(8*x)*exp(-x) axis x1y2"))
output)
(png-from-file
(placing-tics-on-second-y-axis-plot
"images/placing-tics-on-second-y-axis-plot.png"))
Out[38]:
In [39]:
(defun adjusting-the-tic-size-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output
:terminal '(pngcairo)) ;
(gp :set :tics '(scale 3))
(gp :set :mxtics 4)
(gp :set :mytics 4)
(func-plot "[0:4*pi] sin(x)/x"))
output)
(png-from-file
(adjusting-the-tic-size-plot "images/adjusting-the-tic-size-plot.png"))
Out[39]:
In [40]:
(defun removing-all-tics-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :unset :tics)
(func-plot "[0:2*pi] x**2*sin(x)"))
output)
(png-from-file (removing-all-tics-plot "images/removing-all-tics-plot.png"))
Out[40]:
In [41]:
(defun defining-tic-values-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output
:terminal '(pngcairo)) ;
(gp :unset :border)
(gp :set :xtics) '(axis nomirror)
(gp :set :ytics '(axis nomirror))
(gp :set :zeroaxis)
(gp :set :polar)
(gp :set :samples 500)
(gp :set :grid)
(gp :unset :raxis)
(gp :set :xtics '(15 #\, 2 #\, 30))
(gp :set :ytics) '(10 #\, 2 #\, 20)
(func-plot "[0:12*pi] t"))
output)
(png-from-file
(defining-tic-values-plot "images/defining-tic-values-plot.png"))
Out[41]:
In [42]:
(defun making-the-tics-stick-out-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :set :tics 'out)
(gp :set :grid 'nil)
(func-plot "[-3:3] tanh(x)"))
output)
(png-from-file
(making-the-tics-stick-out-plot "images/making-the-tics-stick-out-plot.png"))
Out[42]:
In [43]:
(defun setting-manual-tics-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :set :xtics
'(#\( "pi" pi #\, "pi/2" pi/2 #\, "2pi" 2*pi #\, "3pi/2" 3*pi/2 #\, "0"
0 #\)))
(func-plot "[0:2*pi] sin(x)"))
output)
(png-from-file
(setting-manual-tics-plot "images/setting-manual-tics-plot.png"))
Out[43]:
In [44]:
(defun setting-manual-tics-enhanced-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo enhanced))
(gp :set :xtics
'(#\( "{/symbol p}" pi #\, "{/symbol p}/2" pi/2 #\, "2{/symbol p}" 2*pi
#\, "3{/symbol p}/2" 3*pi/2 #\, "0" 0 #\)))
(func-plot "[0:2*pi] sin(x)"))
output)
(png-from-file
(setting-manual-tics-enhanced-plot
"images/setting-manual-tics-enhanced-plot.png")
)
Out[44]:
In [45]:
(defun setting-manual-tics-enhanced-add-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo enhanced))
(gp :set :xtics '(add #\( "{/symbol p}" pi #\, "2{/symbol p}" 2*pi #\)))
(func-plot "[0:2*pi] sin(x)"))
output)
(png-from-file
(setting-manual-tics-enhanced-add-plot
"images/setting-manual-tics-enhanced-add-plot.png"))
Out[45]:
In [46]:
(defun plotting-with-dates-and-times-plot (output)
(let ((s
(clml.utility.data:fetch
"https://mmaul.github.io/clml.data/sample/traffic-balance.csv"
:stream t)))
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output)
(gp :set :title "interface utilization")
(gp :set :xdata 'time)
(gp :set :timefmt "%y/%m/%d:%h:%m:%s")
(gp :set :format '(x "%d/%m/%y"))
(gp :set :xtics '(rotate))
(gp :unset :key)
(plot
(lambda ()
(loop for row across (subseq
(clml.utility.csv:read-csv-stream s :type-spec
'(string float
float
float
float
float
float))
0 1000)
do (format t "~&~a ~{~a~^ ~}"
(map 'string
(lambda (x)
(if (equal #\ x)
#\:
x))
(elt row 0))
(coerce (subseq row 1) 'list))))
:using '(0 5) :with 'lines)))
output)
(png-from-file
(plotting-with-dates-and-times-plot
"images/plotting-with-dates-and-times-plot.png"))
Out[46]:
The plot below requires gnuplot version 5
In [47]:
(defun changing-language-used-in-labels-plot (output)
(let ((s
(clml.utility.data:fetch
"https://mmaul.github.io/clml.data/sample/traffic-balance.csv"
:stream t)))
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output)
(gp :set :title "interface utilization")
(gp :set :xdata 'time)
(gp :set :locale "es_ES.utf8")
(gp :set :timefmt "%y/%m/%d:%h:%m:%s")
(gp :set :format '(x "%b %y"))
(gp :set :xtics '(rotate))
(gp :unset :key)
(plot
(lambda ()
(loop for row across (subseq
(clml.utility.csv:read-csv-stream s :type-spec
'(string float
float
float
float
float
float))
0 10)
do (format t "~&~a ~{~a~^ ~}"
(map 'string
(lambda (x)
(if (equal #\ x)
#\:
x))
(elt row 0))
(coerce (subseq row 1) 'list))))
:using '(0 5) :with 'lines)))
output)
(png-from-file
(changing-language-used-in-labels-plot
"images/changing-language-used-in-labels-plot.png"))
Out[47]:
The plot below requires gnuplot version 5
In [48]:
(defun using-european-style-decimals-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :set :decimalsign '(locale "es_ES.utf8"))
(func-plot "[0:1] x**3"))
output)
(png-from-file
(using-european-style-decimals-plot
"images/using-european-style-decimals-plot.png"))
Out[48]:
In [49]:
(defun formatting-tic-labels-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :set :ylabel "output power")
(gp :set :format '(y "p = %.2f watts"))
(gp :set :format '(x "%.3f%%"))
(gp :set :rmargin 6)
(func-plot "[0:1] x**3"))
output)
(png-from-file
(formatting-tic-labels-plot "images/formatting-tic-labels-plot.png"))
Out[49]:
In [50]:
(defun tic-no-labels-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo enhanced))
(gp :set :ylabel "output power")
(gp :set :format '(y ""))
(gp :set :format '(x ""))
(gp :set :grid)
(gp :unset :key)
(func-plot "x, -x"))
output)
(png-from-file (tic-no-labels-plot "images/tic-no-labels-plot.png"))
Out[50]:
5 Combining Multiple Plots¶
In [51]:
(defun multiplot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title "multiplot"
:multiplot (list "layout 2,2 columnsfirst title \"multiplot\""))
(func-plot "besj0(x)")
(func-plot "besj1(x)")
(func-plot "besy0(x)")
(func-plot "besy1(x)")
(gp :unset 'multiplot))
output)
(png-from-file (multiplot "images/multiplot.png"))
Out[51]:
In [52]:
(defun stem-multiplot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :title "candlestick"
:xrange (list "[0:5]") :yrange (list "[0:25]") :multiplot
(list "layout 2,1 columnsfirst title \"multiplot stem\""))
(plot
(lambda ()
(loop for r in '((1 5 7 10 15 24 0.3 "'quick'")
(2 6 8 11 16 23 0.4 "'fox'")
(3 5 7 11 17 22 0.5 "'lazy'")
(4 6 9 10 18 21 0.3 "'dog'"))
do (format t "~&~{~^~a ~}" r)))
:using '(1 6 "xticlabels(8)") :with '(impulses lw 4 title "counts") :using
'(1 6) :with '(points pt 9 notitle))
(plot
(lambda ()
(loop for r in '((1 5 7 10 15 24 0.3 "'quick'")
(2 6 8 11 16 23 0.4 "'fox'")
(3 5 7 11 17 22 0.5 "'lazy'")
(4 6 9 10 18 21 0.3 "'dog'"))
do (format t "~&~{~^~a ~}" r)))
:using '(1 6 "xticlabels(8)") :with '(impulses lw 4 title "counts") :using
'(1 6) :with '(points pt 9 notitle))
(gp :unset 'multiplot))
output)
(png-from-file (stem-multiplot "images/stem-multiplot.png"))
Out[52]:
In [53]:
(defun plot-positioning-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :multiplot (list ""))
(gp :set :title "o o")
(gp :set :polar 'nil)
(gp :set :size (list "1,.5"))
(gp :set :border 0)
(gp :unset :key)
(gp :unset :tics)
(gp :unset :raxis)
(func-plot "[pi:2*pi] -1" :lw 5)
(gp :set :title "###################")
(gp :set :origin (list "0,.5"))
(gp :set :size (list ".5, .5"))
(func-plot "-2*pi" :lw (list "2, .2") :with 'filledcurves)
(gp :set :origin (list ".5, .5"))
(gp :set :title "###################")
(func-plot "1" :lw (list "2, .2") :with 'filledcurves)
(gp :unset 'multiplot))
output)
(png-from-file (plot-positioning-plot "images/plot-positioning.png"))
Out[53]:
In [54]:
(defun inset-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(:pngcairo) :multiplot (list ""))
(gp :set :object '(ellipse center 0.13 #\, 0 size 0.4 #\, 4))
(gp :set :arrow '(from 0.1 #\, 2.1 to screen 0.22 #\, 0.4 front lt 3))
(gp :set :samples 1000)
(gp :set :grid 'nil)
(gp :set :xtics 0.4)
(gp :set :ytics 4)
(func-plot "[0:2*pi] exp(x)*sin(1/x)")
(gp :set :origin '(0.2 #\, 0.4))
(gp :set :size '(0.25 #\, 0.25))
(gp :clear)
(gp :unset :grid)
(gp :unset :object)
(gp :unset :arrow)
(gp :set :xtics 0.1)
(gp :set :ytics 0.5)
(gp :set :bmargin 1)
(gp :set :tmargin 1)
(gp :set :lmargin 3)
(gp :set :rmargin 1)
(func-plot "[0:.2] exp(x)*sin(1/x)")
(gp :unset 'multiplot))
output)
(png-from-file (inset-plot "images/inset-plot.png"))
Out[54]:
In [55]:
(defun multiplot-w-labels-and-arrows (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo) :multiplot
'(layout 2 #\, 2 title "derivatives of sin(x)" font
"times-roman, 22"))
(gp :unset :key)
(gp :set :xrange '("[-pi:pi]"))
(gp :set :arrow '(from 0.1 #\, 2.1 to screen 0.22 #\, 0.4 front lt 3))
(gp :set :style
'(arrow 1 head filled size screen 0.03 #\, 15 #\, 135 lt 2 lw 2))
(gp :set :arrow
'(1 from screen 0.45 #\, 0.84 to screen 0.65 #\, 0.84 arrowstyle 1))
(gp :set :arrow
'(2 from screen 0.87 #\, 0.64 to screen 0.87 #\, 0.3 arrowstyle 1))
(gp :set :arrow
'(3 from screen 0.7 #\, 0.15 to screen 0.4 #\, 0.15 arrowstyle 1))
(gp :set :arrow
'(4 from screen 0.35 #\, 0.35 to screen 0.35 #\, 0.7 arrowstyle 1))
(gp :set :title "sin(x)")
(func-plot "sin(x)")
(gp :set :title "sin'(x) = cos(x)")
(func-plot "cos(x)")
(gp :set :title "sin'''(x) = cos''(x) = -sin'(x) = -cos(x)")
(func-plot "-cos(x)")
(gp :set :title "sin''(x) = cos'(x) = -sin(x)")
(func-plot "-sin(x)")
(gp :unset 'multiplot))
output)
(png-from-file
(multiplot-w-labels-and-arrows "images/multiplot-w-labels-and-arrows.png"))
Out[55]:
7 Proramming gnuplot and dealing with data¶
In [56]:
(defun plotting-on-sub-intervals-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output)
(gp :set :samples 2000)
(gp "f(x) = x < 0 ? sin(x) : NaN")
(gp "g(x) = x >= 0 ? exp(-x/5.)*sin(x) : NaN")
(func-plot "[-20 : 20] f(x), g(x)")
)
output)
(png-from-file
(plotting-on-sub-intervals-plot "images/plotting-on-sub-intervals-plot.png"))
Out[56]:
Not working needs to repeat data section again for smoothing, but look for fix in the future
In [57]:
(defun smoothing-your-data-plot (output)
(let ((s
(clml.utility.data:fetch
"http://www.google.com/finance/getprices?i=60&p=10d&f=d,o,h,l,c,v&df=cpct&q=ibm"
:stream t)))
(clml.utility.data:process-finance-header s)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output)
(plot
(lambda ()
(loop for row across (subseq
(clml.utility.csv:read-csv-stream s :type-spec
'(float float
float
float
float
float))
0 300)
do (format t "~&~{~a~^ ~}" (coerce row 'list))))
:using '(0 4) :with
'(lines lw 0.5 notitle #\, ("''") smooth bezier lw 4 title
"bezier smoothed"))))
output)
;(png-from-file
; (smoothing-your-data-plot "images/smoothing-your-data-plot.png"))
Out[57]:
Not working needs to repeat data section again for smoothing, but look for fix in the future
In [58]:
(defun fit-functions-to-your-data-plot (output)
(let ((s
(clml.utility.data:fetch
"http://www.google.com/finance/getprices?i=60&p=10d&f=d,o,h,l,c,v&df=cpct&q=ibm"
:stream t)))
(clml.utility.data:process-finance-header s)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output)
(format t "~&f(x) = a*sin(b*x)~%")
(format t "~&fit f(x) 'rs.dat' via a, b~%")
(plot
(lambda ()
(loop for row across (subseq
(clml.utility.csv:read-csv-stream s :type-spec
'(float float
float
float
float
float))
0 300)
do (format t "~&~{~a~^ ~}" (coerce row 'list))))
:using '(0 4) :with
'(lines lw 0.5 notitle #\, f (x) lw 4 title " fit by gnuplot"))))
output)
;(png-from-file
; (fit-functions-to-your-data-plot
; "images/fit-functions-to-your-data-plot.png"))
Out[58]:
Not working needs to repeat data section again for smoothing, but look for fix in the future
In [59]:
(defun using-kdensity-to-improve-on-histograms-plot (output)
(let ((s
(clml.utility.data:fetch
"http://www.google.com/finance/getprices?i=60&p=10d&f=d,o,h,l,c,v&df=cpct&q=ibm"
:stream t)))
(clml.utility.data:process-finance-header s)
(with-plots (*standard-output* :debug nil)
(gp-setup :terminal '(pngcairo) :output output)
(format t "~&f(x) = a*sin(b*x)~%")
(format t "~&fit f(x) 'rs.dat' via a, b~%")
(plot
(lambda ()
(loop for row across (subseq
(clml.utility.csv:read-csv-stream s :type-spec
'(float float
float
float
float
float))
0 300)
do (format t "~&~{~a~^ ~}" (coerce row 'list))))
:using '(("1:(.001)") smooth kdensity))))
output)
;(png-from-file
; (using-kdensity-to-improve-on-histograms-plot
; "images/using-kdensity-to-improve-on-histograms-plot.png"))
Out[59]:
8 Third dimension¶
In [60]:
(defun 3d-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :unset :key)
(gp :set :isosamples 40)
(gp :set :title "j_0(r^2)")
(gp :set :xrange '("[-4:4]"))
(gp :set :yrange '("[-4:4]"))
(gp :set :ztics 1)
(gp :set :view '(29 #\, 53))
(func-splot "besj0(x**2+y**2)")
(format t "~&replot~%"))
output)
(png-from-file (3d-plot "images/3d-plot.png"))
Out[60]:
In [61]:
(defun pm3d-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :unset :key)
(gp :set :isosamples 40)
(gp :set :title "j_0(r^2)")
(gp :set :xrange '("[-4:4]"))
(gp :set :yrange '("[-4:4]"))
(gp :set :ztics 1)
(gp :set :view '(29 #\, 53))
(gp :set :pm3d 'nil)
(gp :unset :surface)
(func-splot "besj0(x**2+y**2)")
(format t "~&replot~%"))
output)
(png-from-file (pm3d-plot "images/pm3d-plot.png"))
Out[61]:
In [62]:
(defun pm3d2-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :unset :key)
(gp :set :isosamples 40)
(gp :set :title "j_0(r^2)")
(gp :set :xrange '("[-4:4]"))
(gp :set :yrange '("[-4:4]"))
(gp :set :ztics 1)
(gp :unset :surf)
(gp :set :view '(29 #\, 53))
(gp :set :style '(line 1 lt 4 lw 0.5))
(gp :set :pm3d '(at s hidden3d 1))
(func-splot "besj0(x**2+y**2)")
(format t "~&replot~%"))
output)
(png-from-file (pm3d2-plot "images/pm3d2-plot.png"))
Out[62]:
In [63]:
(defun contour-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :unset :key)
(gp :set :cntrparam '(levels 10))
(gp :set :contour 'base)
(gp :unset :sur)
(gp :set :view 'map)
(gp :set :xrange '("[-4:4]"))
(gp :set :yrange '("[-4:4]"))
(gp :set :iso 100)
(gp :set :samp 100)
(gp :set :key 'rmargin)
(func-splot "sin(x)+cos(2*y)"))
output)
(png-from-file (contour-plot "images/contour-plot.png"))
Out[63]:
In [64]:
(defun contour-both (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :unset :key)
(gp :set :cntrparam '(levels 10))
(gp :set :style '(line 1 linecolor rgb "'#cccccc'"))
(gp :set :contour 'both)
(gp :set :hidd 'nil)
(gp :set :iso 50)
(gp :set :samp 50)
(gp :set :yrange '("[-4:4]"))
(gp :set :xrange '("[-4:4]"))
(func-splot "sin(x)+cos(2*y)" :with '(lines linestyle 1)))
output)
(png-from-file (contour-both "images/contour-both.png"))
Out[64]:
In [65]:
(defun vector-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :unset :key)
(gp :set :xrange '("[0:pi]"))
(gp :set :yrange '("[0:pi]"))
(gp :set :iso 20)
(gp :set :samp 20)
(gp :a := .2)
(gp :set :iso 50)
(gp :set :samp 50)
(gp :set :yrange '("[-4:4]"))
(gp :set :xrange '("[-4:4]"))
(func-plot "'++'" :using '(1 2 "(-a*sin($1)*cos($2)):(a*cos($1)*sin($2))")
:with '(vec size 0.06 #\, 15 filled)))
output)
(png-from-file (vector-plot "images/vector-plot.png"))
Out[65]:
In [66]:
(defun vector-contour-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :unset :key)
(gp :set :xrange '("[-pi:pi]"))
(gp :set :yrange '("[-pi:pi]"))
(gp :set :zrange '("[-pi:pi]"))
(gp :set :iso 20)
(gp :set :samp 20)
(gp :set :ztics 1.5)
(gp :set :view '(37 #\, 300))
(gp :a := .9)
(func-splot "'++'" :using
(list 1 2 "(2*a*cos($2)*sin($1))" "(-a*sin($1)*cos($2))"
"(a*cos($1)*sin($2))" "(a*cos($1))")
:with '(vec)))
output)
(png-from-file (vector-contour-plot "images/vector-contour-plot.png"))
Out[66]:
In [67]:
(defun image-map-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :unset :key)
(gp :set :xrange '("[-4:4]"))
(gp :set :yrange '("[-4:4]"))
(gp :set :iso 100)
(gp :set :samp 100)
(gp :set :view 'map)
(gp :set :pm3d '(at b))
(gp :set :palette '("defined (0 'black', 1 'gold')"))
(gp :set :logscale 'cb)
(gp :unset :surface)
(func-splot "sin(y**2+x**2) - cos(x**2)"))
output)
(png-from-file (image-map-plot "images/image-map-plot.png"))
Out[67]:
In [68]:
(defun contours-image-map-plot(output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output
:terminal '(pngcairo))
(gp :unset :key)
(gp :unset :surface)
(gp :set :xrange '("[0:pi]"))
(gp :set :yrange '("[0:pi]"))
(gp :set :iso 100)
(gp :set :samp 100)
(gp :set :cntrparam '(levels 10))
(gp :set :view 'map)
(gp :set :pm3d '(at b))
(gp :set :contour 'base)
;:palette '("defined (0 'black', 1 'gold')")
;:logscale 'cb
;:cbrange "[a : b]"
;(gp :unset :surface )
(func-splot "'++'"
:using '(1 2 "($1**2-$2**2)" "(sin($1**2+$2**2))")
:with '( lines lw 2)))
output)
(png-from-file (contours-image-map-plot "images/contours-image-map-plot.png"))
Out[68]:
In [69]:
(defun surface-w-images-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :unset :key)
(gp :set :xrange '("[-pi:pi]"))
(gp :set :yrange '("[-pi:pi]"))
(gp :set :iso 40)
(gp :set :samp 40)
(gp :set :hidden 'front)
(format t "~&f(x,y) = sin(x)*cos(y)~%")
(func-splot "f(x,y)" :with (list "pm3d at b, f(x,y)") :with 'lines))
output)
(png-from-file (surface-w-images-plot "images/surface-w-images-plot.png"))
Out[69]:
In [70]:
(defun path-in-3d-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :set :xtics 0.4)
(gp :set :ytics 0.4)
(gp :set :parametric 'nil)
(gp :set :samp 40)
(gp :set :ztics 1)
(func-splot "cos(u),sin(3*u),cos(5*u) lw 2"))
output)
(png-from-file (path-in-3d-plot "images/path-in-3d-plot.png"))
Out[70]:
In [71]:
(defun parametric-surfaces (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :set :param 'nil)
(gp :set :iso 50)
(gp :set :ztics 0.5)
(gp :set :xtics 0.4)
(gp :set :ytics 0.4)
(gp :set :urange '("[-pi:pi]"))
(gp :set :vrange '("[-pi:pi]"))
(gp :set :hidd 'nil)
(func-splot "cos(u)*cos(v), sin(u)*cos(v), sin(u)"))
output)
(png-from-file (parametric-surfaces "images/parametric-surfaces.png"))
Out[71]:
In [72]:
(defun parametric-surfaces-depthorder (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :set :param 'nil)
(gp :set :iso 50)
(gp :set :ztics 0.5)
(gp :set :xtics 0.4)
(gp :set :ytics 0.4)
(gp :set :urange '("[-pi:pi]"))
(gp :set :vrange '("[-pi:pi]"))
(gp :set :hidd 'nil)
(gp :set :pm3d 'depthorder)
(func-splot "cos(u)*cos(v), sin(u)*cos(v), sin(u)"))
output)
(png-from-file
(parametric-surfaces-depthorder "images/parametric-surfaces-depthorder.png"))
Out[72]:
10 Surveying Special Topics¶
In [73]:
(defun rotate-labels-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo))
(gp :set :xrange '("[0:2*pi]"))
(gp :set :xtics '(0 #\, ("0.12*pi") #\, ("2*pi")))
(gp :set :xtics '(format "%.10f"))
(gp :set :xtics '(rotate by -30))
(gp :set :rmargin 10)
(func-plot "sin(x)"))
output)
(png-from-file (rotate-labels-plot "images/rotate-labels-plot.png"))
Out[73]:
In [74]:
(defun labels-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo enhanced size 750 #\, 1000)
:title "label plot" :xrange (list "[0:15]") :yrange
(list "[6:10]"))
(gp :unset :key)
(gp :set :rmargin 10)
(gp :set :xlabel '("'one thing'"))
(gp :set :ylabel '("'another'"))
(format t
"~&countryname(string,size) = sprintf(\"{/=%d %s}\", size, string)~%")
(plot
(lambda ()
(loop for r in '((1 5 7 10 15 24 0.3 "'quick'")
(2 10 8 11 20 23 0.4 "'fox'")
(3 9 7 11 7 22 0.5 "'lazy'")
(4 6 9 10 40 21 0.3 "'dog'"))
do (format t "~&~{~^~a ~}" r)))
:using '(2 3 ("(countryname(stringcolumn(8),$5))")) :with 'labels))
output)
(png-from-file (labels-plot "images/labels-plot.png"))
Out[74]:
In [75]:
(defun labeled-contour-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output
:terminal '(pngcairo)
:view :map
:contour :base
:xrange '("[0:5]")
:yrange '("[0:5]")
:style '( fill solid)
:for '("[n = 1:4]" cntrparam levels discrete ("n**2"))
:for '("[n = 1:4]" object n circle at n #\, n size 0.2 front fillcolor rgb ("'#ffffff'") lw 0)
:for '("[n = 1:4]" label n
("sprintf(\"%d\", n**2)") at n #\, n
center front)
)
(gp :unset :key)
(gp :unset :surf)
(gp :unset :clabel)
(func-splot "x*y")
)output)
(png-from-file (labeled-contour-plot "images/labeled-contour-plot.png"))
Out[75]:
In [76]:
(defun softened-axis-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo) :style
'(line 2 lc rgb ("'#999999'") lt 0 ("#grid")) :style
'(line 1 lc rgb ("'#999999'") lt 1 ("#border")) :grid
'(linestyle 2) :border '(linestyle 1))
(gp :unset :key)
(gp :unset :surf)
(gp :unset :clabel)
(func-plot "sin(x)/x" :lw 2))
output)
(png-from-file (softened-axis-plot "images/softened-axis-plot.png"))
Out[76]:
In [77]:
(defun arrow-axis-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo) :arrow
'(from graph 0 #\, 1 to graph 0 #\, 1.1 filled) :arrow
'(from graph 1 #\, 0 to graph 1.1 #\, 0 filled) :tmargin 5
:rmargin 20 :border 3 :tics 'nomirror :grid 'nil)
(func-plot "sin(x)/x" :lw 2))
output)
(png-from-file (arrow-axis-plot "images/arrow-axis-plot.png"))
Out[77]:
In [78]:
(defun plotting-with-images-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo) :title
'("\"elbonian warhead inventory\"" font "times,24" offset 0 #\,
2)
:xlabel '("" "year" font "times,12") :xlabel
'("" "year" font "times,16") :ylabel
'("" "number of operational warheads\\n" font "times,14") :tics
'nomirror :xtics '(scale 0) :for
'("[n = 0 : 3]" ytics ("(sprintf(\"%d\" , (n*10)) 200*n)")) :for
'("[n = 0 : 3]" xtics ("(sprintf(\"%d\" , 1990+10*n) 91+200*n)"))
:xtics '(add ("(\"2020\\n(projected)\" 691)")) :cbrange
'("[0:250]") :pal '(gray positive) :view 'map)
(gp :unset :key)
(gp :unset :colorbox)
(func-splot
" for [n = 0 : 3] 'images/bomb.png' binary filetype=auto center=(91+200*n, 300*0.2*(1+n), 0) dy=0.2*(1+n) with image"))
output)
(png-from-file
(plotting-with-images-plot "images/plotting-with-images-plot.png"))
Out[78]:
In [79]:
(defun broken-axis-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo font "times,8") :border 2
:zeroaxis '(lt -1) :ytics '(axis nomirror) :xtics
'(axis nomirror) :tics 'out :samp 5000 :xrange '("[0.007:.1]")
:yrange '("[-.05:1]") :xtics '(0.005 #\, 0.005 #\, 0.025)
;:object '("1 fillstyle pattern 5 lw 0")
:key '(font "times,12" samplen -1))
(gp :unset :key)
(gp :unset :surf)
(gp :unset :clabel)
(gp "s" := 0.2)
(gp "b" := 0.03)
(gp "w" := 0.001)
(gp "h" := 1)
(gp "h(x)" := '("sin(1/x)**2"))
(gp "f(x)" := '("(x<b?h(x):NaN)"))
(gp "g(x)" := '("(x>b?h(x+s):NaN)"))
(gp :set :for
'("[n = 4 : 10] xtics add(sprintf(\"%.3f\", s + n/100.0) n/100.0)"))
;(gp :set :object '("1 polygon from first b-w,-h to b+w,-h to b+w,h to b-w,h to b-w,-h front"))
(func-plot " f(x) title \"sin(1/x)\" lt -1, g(x) notitle lt -1"))
output)
(png-from-file (broken-axis-plot "images/broken-axis-plot.png"))
Out[79]:
In [80]:
(defun fit-grid-to-data-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo font "times,8") :samples 6
:key '(top left) :for
'("[n = 1 : 4] arrow from first n, 0 to first n, n**2 back nohead lt 7")
:for
'("[n = 1 : 4] arrow from first 0, n**2 to first n, n**2 back nohead lt 7")
:for '("[n = 0 : 5] ytics (n**2)"))
(func-plot "[0:5] x**2 with linespoints pt 7 ps 3"))
output)
(png-from-file (fit-grid-to-data-plot "images/fit-grid-to-data-plot.png"))
Out[80]:
In [81]:
(defun coloring-axis-plot (output)
(with-plots (*standard-output* :debug nil)
(gp-setup :output output :terminal '(pngcairo font "times,8") :multiplot
(list "") :tmargin '(at screen 0.8) :rmargin '(at screen 0.8)
:lmargin '(at screen 0.2) :bmargin '(at screen 0.2) :tics
'nomirror :xrange '("[0:5]") :border 8 :border '(lt 4) :y2tics
'(0 #\, 25 #\, 150) :y2label '("\"exp(x)\"" textcolor lt 4))
(gp :unset :key)
(gp :unset :ytics)
(func-plot "exp(x) axis x1y2 lt 4")
(gp :set :border 7)
(gp :set :border '(lt -1))
(gp :set :ylabel '("\"1/x\"" textcolor lt -1))
(gp :unset :y2label)
(gp :set :ytics '(0 #\, 1 #\, 5))
(gp :set :tics 'nomirror)
(gp :unset :y2tics)
(gp :set :yrange '("[0:5]"))
(gp :unset :y2tics)
(func-plot "1/x" :lt -1)
(gp :unset :multiplot))
output)
(png-from-file (coloring-axis-plot "images/coloring-axis-plot.png"))
Out[81]:
In [ ]: